Drift tube borehole muon detector system, apparatus, and method for muon radiography and tomography

ABSTRACT

A borehole muon detector for muon radiography or geotomography is provided, the borehole muon detector including a substantially cylindrical housing, which defines a bore, a pair of end caps, each end cap sealing an end of the cylindrical housing and a plurality of sealed drift tubes which are longitudinally disposed in the bore of the housing to form a bundle of drift tubes, wherein each sealed drift tube comprises: a centrally located anode wire disposed on a longitudinal axis; an inner surface which is coated with a cathode coating, the cathode coating divided into a first cathode pad and a second cathode pad by a Vernier pattern; and a timer in electrical communication with the anode wire for measuring a drift time. A system and a method are also provided.

REFERENCE TO RELATED APPLICATIONS

This application is a continuation of U.S. application Ser. No.16/975,499 having a filing (371) date of 25 Aug. 2020, which is in turna 371 national phase application corresponding to Patent CooperationTreaty (PCT) application No. PCT/CA2020/000036 having an internationalfiling date of 20 Mar. 2020, which is turn claims priority to and thebenefit under 35 USC 119 in relation to U.S. application Ser. No.62/821,448 filed 21 Mar. 2019. All of the applications referred to inthis paragraph are hereby incorporated herein by reference.

FIELD

The present technology relates generally to a system, apparatus, andmethod for tracking cosmic ray muons through an underground sensor inorder to develop an image of subsurface density above the sensor (muonradiography), and to use multiple sensors to build a 3D model of density(muon tomography).

BACKGROUND

Exploitation of underground resources (e.g., mineral deposits, oilreservoirs) employs varied geophysical methods to detect, image, andmonitor underground regions of interest. One such method is based on atechnique known as muon radiography. Muons are elementary particlesproduced in high energy nuclear interactions that are initiated bycosmic rays in the upper atmosphere. The muons travel at nearly thespeed of light and, depending on their energy, can penetrate deep intothe earth. The rate at which the muons lose energy in matter isdependent on the properties of the medium, and in particular on thedensity of matter.

If the rate of muons at all energies and angles at the surface of theearth is known, it can be used to predict the rate of muons passingthrough a sensor underground, assuming an underlying geological model.This relationship can be inverted to develop a geological model based onthe observed rate of muons underground. Muon sensors record thetrajectory of muons that pass through them. Over a period of exposuretime, the rate of muons passing from any given direction can bemeasured. From this observed rate, one can infer the average densityalong a ray path from the sensor towards the surface.

Thus, from a single muon sensor, a map, or image, of average subsurfacedensity above the sensor can be ascertained.

The resolution of the density map is determined by the angularresolution of the muon sensor and by the lateral scattering of muons asthey pass through matter. The lateral scattering of muons through fiftymeters of rock is about twenty-five milliradians, on average. Theaverage scattering angle decreases with increasing depth because themean energy of muons that survive increases, and high energy muons aredeflected less than low energy ones.

Density maps from multiple sensor locations (either using the samesensor exposed to muon flux in multiple locations, or else from multiplesensors) can be combined to construct a 3D model of the undergrounddensity distribution. This is known as muon tomography. This process isachieved by applying inversion algorithms to the set of sensor images,and proceeds analogously to computed tomography for X-ray scanners. Themore varied locations from which muon rates can be measured, the morestrongly a corresponding 3D density model of the subsurface above thesensors can be constrained.

Anomalies in the subsurface density distribution are often indicative ofvaluable resources such as dense mineral deposits (e.g., volcanogenicmassive sulfide deposits, uranium deposits), oil reservoirs, or othergeological markers commonly associated with such resources (e.g.,faults). Such anomalies can be identified either in the sensor images(radiography), or with improved spatial information in the 3Ddistributions (tomography). Identification of density anomalies is alsorelevant to geotechnical applications such as infrastructure monitoring.

United States Patent Application Publication No. 20200018863 discloses amuon detector system capable of determining muon direction and flighttrajectory or path is disclosed. The muon detector system includesscintillators for determining muon direction, and an array of muondetectors arranged in orthogonal layers of drift tubes for determiningflight trajectory. The system can be used for tomographic and telescopicmode imaging and may be used for imaging concealed and/or subterraneanobjects. It is not, however, suitable for borehole applications.

United States Patent Application Publication No. 20140319365 disclosesmethods, systems, and devices are disclosed for inspecting materials ina vehicle or object. In one aspect, a system for muon tomographydetection includes a first and second housing structure each including afirst array and second array of muon detection sensors, respectively,the first housing structure positioned opposite the second at a fixedheight to form a detection region to contain a target object, in whichthe muon detection sensors measure positions and directions of muonspassing through the first array to the detection region and passing fromthe detection region through the second array; support structures toposition the first housing structure at the fixed height; and aprocessing unit to receive data from the muon detection sensors andanalyze scattering behaviors of the muons in materials of the targetobject to obtain a tomographic profile or spatial distribution ofscattering centers within the detection region. In one exemplaryembodiment, each of the upper and lower detection units includes threesets of horizontally X direction-oriented detector arrays alternativelyvertically interspersed with three sets of horizontally Y directionoriented arrays, with the X and Y directions being perpendicular to eachother. The size and shape of the system precludes it from being usefulfor borehole applications. It is also not applicable for anything largerthan a few meters thick because (a) the system would need to beunmanageably large and (b) the solid angle coverage gets smaller andsmaller as the separation between the trackers increases.

United States Patent Application Publication No. 20080128604 discloses asystem and method configured for detecting cosmic ray muon (CRM) fluxalong a variety of trajectories through a subterranean test region,collecting the muon detection data, and processing the data to form athree-dimensional density distribution image corresponding to the testregion. In at least one embodiment, the system may be used foridentifying concentrations of higher (or lower) density mineral depositsor other geological structures or formations well below the earth's (orocean floor) surface. The system may be utilized for imaging geologicalmaterials and structures of higher and/or lower density in a test regionhaving a depth of several kilometers or more. The time projectionchamber has a restricted length because of the requirement that theionization trail drift needs to be precisely controlled, and this isvery difficult for a long cylinder.

WIPO Publication No. WO/2009/002602 discloses techniques, apparatus andsystems for detecting particles such as muons for imaging applications.Subtraction techniques are described to enhance the processing of themuon tomography data. Each set of position sensitive detectors comprisesa first double-layer 220 of drift tubes 204 arranged in the X directionand a second double-layer 221 of drift tubes 204 arranged in the Ydirection. The size and shape of the system precludes it from beinguseful for borehole applications. Orthogonal planar layers cannot beused in a borehole system because the transverse dimension is limited tobe less than 10 cm in a borehole system.

What is needed is a muon detector that is suitable for boreholeapplications. It would be preferable if it had a bundle of drift tubes.It would be further preferable if the borehole muon detector did notrequire orthogonally disposed drift tubes and instead, had a Vernierpattern on the cathode of each drift tube to allow precise determinationof the longitudinal coordinate of a muon. It would be further preferableif it included longitudinally disposed scintillators to provide areference time of when a muon impinged upon the detector. It would bepreferable if the borehole muon detector was tubular. It would bepreferable if it was a sealed unit that could withstand high pressureand would be waterproof. It would be most preferable if it could measurethe trajectory of muons in three dimensions.

SUMMARY

The present technology is a borehole muon detector. It has drift tubesthat are longitudinally disposed and in a bundle. Each drift tubeincludes a Vernier pattern on the cathode to allow precise determinationof the longitudinal coordinate of a muon. The borehole muon detectorincludes longitudinally disposed scintillators on an outer surface ofthe bundle of drift tubes to provide a reference time of when a muonimpinged upon the detector. The borehole detector is tubular and is asealed unit that can withstand high pressure. It is waterproof anddustproof. The borehole muon detector can measure the trajectory of amuon in three dimensions.

In one embodiment borehole muon detector for muon radiography orgeotomography is provided, the borehole muon detector including asubstantially cylindrical housing, which defines a bore, a pair of endcaps, each end cap sealing an end of the cylindrical housing and aplurality of sealed drift tubes which are longitudinally disposed in thebore of the housing to form a bundle of drift tubes, wherein each sealeddrift tube comprises: a centrally located anode wire disposed on alongitudinal axis; an inner surface which is coated with a cathodecoating, the cathode coating divided into at least a first cathode padand a second cathode pad by a Vernier pattern; and a printed circuitboard in electrical communication with the anode wire for measuring acurrent.

In the borehole muon detector, there may be at least three drift tubesin the bundle.

The borehole muon detector may further include a plurality ofscintillator members disposed on the bundle.

The borehole muon detector may further include at least one positionsensor.

In the borehole muon detector, the position sensor may be one or more ofan accelerometer, a magnetometer and a gyroscope.

In another embodiment, a system for borehole muon geotomography isprovided, the system including at least one borehole muon detector forplacement in at least one boreholes and a processor in electroniccommunication with each of the borehole muon detectors, each boreholemuon detector including a substantially cylindrical housing, whichdefines a bore, a pair of end caps, each end cap sealing an end of thecylindrical housing and a plurality of sealed drift tubes which arelongitudinally disposed in the bore of the housing to form a bundle ofdrift tubes, wherein each sealed drift tube comprises: a centrallylocated anode wire disposed on a longitudinal axis; an inner surfacewhich is coated with a cathode coating; and printed circuit board inelectrical communication with the anode wire and in electroniccommunication with the processor.

The system may further comprise a plurality of scintillator membersdisposed on at least one bundle.

In the system, each borehole muon detector may further include at leastone position sensor.

In the system, the position sensor may be one or more of anaccelerometer, a magnetometer and a gyroscope.

In the system, the cathode coating may be divided into at least a firstcathode pad and a second cathode pad by a Vernier pattern.

In another embodiment, a method for borehole muon radiography orgeotomography is provided, the method comprising: positioning aplurality of borehole muon detectors in one or more boreholes, eachborehole muon detector including a plurality of sealed drift tubes, eachdrift tube including a centrally located anode wire disposed on alongitudinal axis, an inner surface which is coated with a cathodecoating and timer; generating a radial electric field between the anodewire and the cathode coating; and measuring a plurality of muontrajectories.

In the method, the plurality of muon trajectories may be processed toprovide a two-dimensional image.

In the method, the plurality of muon trajectories may be processed toprovide a plurality of two-dimensional images.

In the method, the plurality of two-dimensional images may be furtherprocessed by applying inversion algorithms to provide at least one threedimensional image.

FIGURES

FIG. 1 is a schematic of a borehole muon detector system of the presenttechnology.

FIG. 2A is a cross sectional view of a borehole muon detector of thesystem of FIG. 1 ; and FIG. 2B is a longitudinal sectional view of aborehole muon detector of the system of FIG. 1 .

FIG. 3 is a perspective longitudinal view of a drift tube of theborehole muon detector of FIGS. 2A and 2B.

FIG. 4 shows an exemplary “Vernier pattern” cathode system on the innerwall of a drift tube.

FIG. 5 is a schematic of the electronics of the borehole muon detectorof FIGS. 2A and 2B.

FIG. 6A is a cross sectional view showing how a muon is tracked througha single drift tube of the borehole muon detector of FIGS. 2A and 2B;FIG. 6B is a longitudinal view showing how a muon is tracked through anumber of drift tubes of the borehole muon detector of FIGS. 2A and 2B.

FIG. 7A is a cross sectional view showing how a muon is tracked throughthe borehole muon detector of FIGS. 2A and 2B; and FIG. 7B is alongitudinal view showing how a muon is tracked through the boreholemuon detector of FIGS. 2A and 2B.

FIGS. 8A and 8B shows sensor images for a simulated cuboid, 100 metersabove a muon sensor situated underground. FIG. 8A shows the image with asliding window filtering algorithm; FIG. 8B shows with image without thesliding window filtering algorithm; and FIGS. 8C and 8D show that eachpixel represents a region of solid angle.

FIG. 9 shows an array of boreholes with borehole muon detectorsinstalled and the detection of a density anomaly using an array ofborehole muon detectors in the borehole muon detector system of FIG. 1 .

DESCRIPTION

Except as otherwise expressly provided, the following rules ofinterpretation apply to this specification (written description andclaims): (a) all words used herein shall be construed to be of suchgender or number (singular or plural) as the circumstances require; (b)the singular terms “a”, “an”, and “the”, as used in the specificationand the appended claims include plural references unless the contextclearly dictates otherwise; (c) the antecedent term “about” applied to arecited range or value denotes an approximation within the deviation inthe range or value known or expected in the art from the measurementsmethod; (d) the words “herein”, “hereby”, “hereof1, “hereto”,“hereinbefore”, and “hereinafter”, and words of similar import, refer tothis specification in its entirety and not to any particular paragraph,claim or other subdivision, unless otherwise specified; (e) descriptiveheadings are for convenience only and shall not control or affect themeaning or construction of any part of the specification; and (f) “or”and “any” are not exclusive and “include” and “including” are notlimiting. Further, the terms “comprising,” “having,” “including,” and“containing” are to be construed as open-ended terms (i.e., meaning“including, but not limited to,”) unless otherwise noted.

Recitation of ranges of values herein are merely intended to serve as ashorthand method of referring individually to each separate valuefalling within the range, unless otherwise indicated herein, and eachseparate value is incorporated into the specification as if it wereindividually recited herein. Where a specific range of values isprovided, it is understood that each intervening value, to the tenth ofthe unit of the lower limit unless the context clearly dictatesotherwise, between the upper and lower limit of that range and any otherstated or intervening value in that stated range, is included therein.All smaller sub ranges are also included. The upper and lower limits ofthese smaller ranges are also included therein, subject to anyspecifically excluded limit in the stated range.

Unless defined otherwise, all technical and scientific terms used hereinhave the same meaning as commonly understood by one of ordinary skill inthe relevant art. Although any methods and materials similar orequivalent to those described herein can also be used, the acceptablemethods and materials are now described.

A borehole muon detector system, generally referred to as 8 is shown inFIG. 1 . It includes one or more borehole muon detectors 10, a computerprocessing unit 12 and a user interface 14. Support cables 16 andelectrical connections 18 extend from the borehole muon detector 10 tothe computer processing unit 12. The borehole muon detector system 8 canbe used in resource exploration and monitoring, and in securityapplications. The system measures the trajectory of muons that passthrough it, from which a differential intensity measurement(radiographic image) can be derived. The differential intensitymeasurements from multiple borehole muon detection systems 8 can beinverted to estimate the three-dimensional distribution of subsurfacedensity (tomographic analysis).

As shown in FIG. 2A, each borehole muon detector 10 includes drift tubes20 which are arranged in a bundle, generally referred to as 22, with aplurality of scintillator members 24 around the periphery 26 of thebundle 22. The drift tubes 20 and the scintillator members 24 areencased in a housing 30, which is tubular. As shown in FIG. 2B, thedrift tubes 20 and the scintillator members 24 are longitudinallydisposed in the bore 32 of the housing 30 and lie overtop the drift tubebundle 22. It can be seen that they are segmented from one another.There are no orthogonally disposed drift tubes 20. The housing 30 has anendcap 34 at each end 36, 38. The endcaps 34 include apertures 40 forthe electrical connections 16 and support cables 18. At one end 36 ofthe housing 30 there is a power input 42 and data output connector 44.The data output connector 44 consists of a ruggedized Ethernet orcoaxial connector. The power input 42 is a suitable ruggedized DCelectrical connector. At the other end 38 the housing 30 there is apower output 46 and data input connector 48. The data input connector 48consists of a ruggedized Ethernet or coaxial connector. The power output46 is a suitable ruggedized DC electrical connector. A support matrix 50locates and retains the drift tubes 20.

As shown in FIG. 3 , each drift tube 20 includes a tubular wall 56 ofrigid material such as plastic, aluminum or carbon fiber, or a Mylar®film that is appropriately glued or supported to maintain rigidity. Thedrift tube 20 has a first end 58 and a second end 60. Each cylindricaldrift tube 20 is a few centimeters in diameter and many centimeters (upto a few meters) in length. The inner surface 62 of each drift tube wall56 has an electrically conductive surface coating 64, which is acathode.

The tubular wall 56 defines a bore 72. The bore 72 is filled with anionizing gas mixture such as argon/carbon dioxide or other appropriatemixture. Centrally located in the bore 72 is an anode wire 74. Thecentral anode wire 74 is an appropriate tungsten/gold/steel wire,affixed in place by plugs 76 at both ends of the tubular wall 60, andmaintained at a high potential of a few thousand volts with respect tothe cathode coating 64. The anode wire 74 has a first end 78 and asecond end 80. It is disposed on the longitudinal axis 82. The drifttube 20 is sealed.

As shown in FIG. 4 , the cathode coating 64 is divided into two separateparts, forming a “Vernier” pattern 66 which extends from the first end36 to the second end 38, to provide a first cathode pad 68 and a secondcathode pad 70. There may be more than two cathode pads 68, 70.

As shown in FIG. 5 , each borehole muon detector 8 includes a printedcircuit board (PCB) 90. The power input 42, data output connector 44,the power output 46 and the data input connector 48 are in electricaland electronic communication with the PCB 90.

The PCB 90 consists of preamplifiers 92, a high voltage generator 94 tosupply potential to the drift tube anode wires 74, application specificintegrated circuits (ASICs) 96, analog to digital converters (ADCs) 98for measuring the voltage on the cathode pads 68, 70 of each drift tube20, the current on the anode wires 74, and time to digital converters(TDCs) 100 for measuring the time of a pulses. Alternatively,“chargedivision” can be used with a resistive anode wire 74 instead of relativetiming on the anode wires 74. The TDCs 100 determine both the drift timeof the ionization and the relative timing along the anode wire 74 of thepulse arrival at each end 36, 38 of the drift tube 20. The latter isused to determine which of the repeated Vernier segments (a pair ofcathode pads 68, 70) is aligned with the ionization event. The PCB 90also includes a field programmable gate array (FPGA) 102, anaccelerometer 104, a magnetometer 106, and a gyroscope 108. The TDCs 100and ADCs 96 are connected to field programmable gate arrays (FPGAs)which process the data from all drift tubes 20 to identify coincidencesconsistent with the passage of a muon, and export this data to an onlinebackend computer 18 for further processing, storage, and pipelining to asurface data collector 1 10. Additional electronics for monitoringtemperature, gas pressure, and other slow controls are embedded withinthe system.

As shown in FIG. 6A, the electric field 120 is perpendicular to thelongitudinal axis 82 and therefore the anode wire 72, as it traversesthe bore 72 from the cathode coating 64 to the anode wire 72. The muons122 create ionization 124 that drifts towards the anode wire 72 on driftpaths 126. The radial position within the drift tube 20 (how close themuon passes by the anode wire 74) is determined by the time it takes forthe closest electrons to drift to the anode wire 74. By measuring thedifference in time or current between the current pulses measured oneach end 78, 80 of the anode wire 74, a rough estimate of the positionalong the anode wire 74 where the muon passed by is achievable. Eachdrift tube 20 measures position along the anode wire 74 (Z measurement),and a drift time relative to a reference clock (may be provided by ascintillator pulse from the scintillator members 24).

As shown in FIG. 6B, the Vernier pattern 66 of the cathode coating 64then allows for a much finer measurement of longitudinal position.Within a single Vernier pattern the fractional mirror charge measured onthe respective cathodes is linearly related to longitudinal position. Ifa muon 122 passes through at least three drift tubes 20 within thehousing 30, the azimuthal angle of the muon trajectory can be accuratelydetermined from the radial measurements. The zenith angle of the muontrajectory 128 with respect to the longitudinal axis 82 of the boreholemuon detector 10 is determined from regression on the longitudinalpositions measured in each of the drift tubes 20. Thus, a borehole muondetector 10 comprised of an array of drift tubes 20 instrumented with aVernier pattern 66 on the cathode coating 64 of the drift tubes 20,allows for precise determination of muons 122 that pass through theborehole muon detector 10. Difference in drift times between multipledrift tubes 20 in a bundle 22 provides X-Y measurements and with X-Y-Zmeasurements from multiple drift tubes 20, the full trajectory can bereconstructed.

FIG. 7A shows how a muon is tracked through the borehole muon detector10 in the transverse plane. A line is fit through all the drift tubes 20in which a pulse is observed, using the isochrones (white circles) toperform a best fit. FIG. 7B shows how a muon is tracked through theborehole muon detector 10 in the longitudinal plane. A line is fit usingthe longitudinal position estimates (white points) coming from timedifferences between first end 58 and second end 60 measurements of eachdrift tube 20 and from the voltage division among the “Vernier pattern”cathodes 64.

FIGS. 8A and 8B show two dimensional radiographic images of muonintensity underneath a simulated 50 m×50 m×20 m cuboid (1 g/cccontrast), 100 meters above a borehole muon detector that is situated400 meters underground. The image is shown with (FIG. 8A) and without(FIG. 8B) a sliding window filtering algorithm applied. Each pixelrepresents a region of solid angle as shown in FIGS. 8C and 8D. The muonintensity shadow shows the location of a density anomaly above theborehole muon detector 10. Three dimensional images can be obtained fromthe two dimensional images by applying inversion algorithms to a set orplurality of two dimensional images.

FIG. 9 is a schematic of a borehole muon detector system 8 with multipleborehole muon detectors 10 arranged in series in each of two boreholes.

In an embodiment, the borehole muon detector system 8 or the boreholemuon detector 10 may determine an occupancy or muon intensity map(radiographic image) of a projected surface above the borehole muondetector system 8 or the borehole muon detector 10 using the muontrajectory for muons that pass through the borehole muon detector system8 or the borehole muon detector 10. The muon intensity in each sectionof solid angle emanating from the borehole muon detector system 8 or theborehole muon detector 10 to the surface is a measure of the averagedensity of the earth within that volume. Regions of lower or higherdensity will have correspondingly higher or lower muon occupancy withinthe corresponding pixels of the radiographic images. The radiographicimages from the borehole muon detector system 8 may facilitatetomographic analysis in order to estimate the three-dimensionaldistribution of subsurface density above the borehole systems.

The invention comprises a number of non-limiting aspects. Non-limitingaspects of the invention provide:

-   -   1. A borehole muon detector for muon radiography or        geotomography, the borehole muon detector including a        substantially cylindrical housing, which defines a bore, a pair        of end caps, each end cap sealing an end of the cylindrical        housing and a plurality of sealed drift tubes which are        longitudinally disposed in the bore of the housing to form a        bundle of drift tubes, wherein each sealed drift tube comprises:        a centrally located anode wire disposed on a longitudinal axis;        an inner surface which is coated with a cathode coating, the        cathode coating divided into at least first cathode pad and a        second cathode pad by a Vernier pattern; and a printed circuit        board in electrical communication with the anode wire for        measuring a current.    -   2. The borehole muon detector of aspect 1, wherein there are at        least three drift tubes in the bundle.    -   3. The borehole muon detector of aspect 2, further including a        plurality of scintillator members disposed on the bundle.    -   4. The borehole muon detector of aspect 2 or 3, further        including at least one position sensor.    -   5. The borehole muon detector of aspect 4, wherein the position        sensor is one or more of an accelerometer, a magnetometer and a        gyroscope.    -   6. A system for borehole muon geotomography, the system        including at least one borehole muon detector for placement in        at least one boreholes and a processor in electronic        communication with each of the borehole muon detectors, each        borehole muon detector including a substantially cylindrical        housing, which defines a bore, a pair of end caps, each end cap        sealing an end of the cylindrical housing and a plurality of        sealed drift tubes which are longitudinally disposed in the bore        of the housing to form a bundle of drift tubes, wherein each        sealed drift tube comprises: a centrally located anode wire        disposed on a longitudinal axis; an inner surface which is        coated with a cathode coating; and printed circuit board in        electrical communication with the anode wire and in electronic        communication with the processor.    -   7. The system of aspect 6, further comprising a plurality of        scintillator members disposed on at least one bundle.    -   8. The system of aspect 6 or 7, wherein each borehole muon        detector further includes at least one position sensor.    -   9. The system of aspect 8, wherein the position sensor is one or        more of an accelerometer, a magnetometer and a gyroscope.    -   10. The system of any one of aspects 6 to 9, wherein the cathode        coating is divided into at least a first cathode pad and a        second cathode pad by a Vernier pattern.    -   11. A method for borehole muon radiography or geotomography, the        method comprising: positioning a plurality of borehole muon        detectors in one or more boreholes, each borehole muon detector        including a plurality of sealed drift tubes, each drift tube        including a centrally located anode wire disposed on a        longitudinal axis, an inner surface which is coated with a        cathode coating and timer; generating a radial electric field        between the anode wire and the cathode coating; and measuring a        plurality of muon trajectories.    -   12. The method of aspect 11, wherein the plurality of muon        trajectories are processed to provide a two dimensional image.    -   13. The method of aspect 11, wherein the plurality of muon        trajectories are processed to provide a plurality of two        dimensional images.    -   14. The method of aspect 13, wherein the plurality of        two-dimensional images are further processed by applying        inversion algorithms to provide at least one three dimensional        image.

While example embodiments have been described in connection with what ispresently considered to be an example of a possible most practicaland/or suitable embodiment, it is to be understood that the descriptionsare not to be limited to the disclosed embodiments, but on the contrary,is intended to cover various modifications and equivalent arrangementsincluded within the spirit and scope of the example embodiment. Thoseskilled in the art will recognize or be able to ascertain using no morethan routine experimentation, many equivalents to the specific exampleembodiments specifically described herein. Such equivalents are intendedto be encompassed in the scope of the claims, if appended hereto orsubsequently filed.

1. A system for borehole muon geotomography, the system including atleast one borehole muon detector for placement in at least one boreholeand a processor in electronic communication with each of the boreholemuon detectors, each borehole muon detector including a substantiallycylindrical housing, which defines a bore, a pair of end caps, each endcap sealing an end of the cylindrical housing and a plurality of sealeddrift tubes longitudinally disposed in the bore of the housing to form abundle of drift tubes, wherein each sealed drift tube comprises: acentrally located anode wire disposed on a longitudinal axis; an innersurface which is coated with a cathode coating; and a printed circuitboard in electrical communication with the anode wire and in electroniccommunication with the processor for determining X-Y measurements in atransverse plane of a passing muon based on drift times of multipledrift tubes in the bundle.
 2. The system of claim 1, wherein determiningthe X-Y measurements in the transverse plane of the passing muon basedon drift times of multiple drift tubes in the bundle comprisesdetermining, for each of the multiple drift tubes, a radial position ofthe muon passing by the anode wire, and, based on the radial positionsfor at least three drift tubes, an azimuthal angle in the transverseplane based on fitting a line through isochrones of the at least threedrift tubes.
 3. The system of claim 1 wherein the processor isconfigured to determine a zenith angle of the muon trajectory withrespect to a longitudinal axis based on regression of longitudinalpositions of the passing muon measured in each of the multiple drifttubes.
 4. The system of claim 1, further comprising a plurality ofscintillator members disposed on at least one bundle.
 5. The system ofclaim 1, wherein each borehole muon detector further includes at leastone position sensor.
 6. The system of claim 5, wherein the positionsensor is one or more of an accelerometer, a magnetometer and agyroscope.
 7. The system of claim 1, wherein the cathode coating isdivided into at least a first cathode pad and a second cathode pad by aVernier pattern.
 8. A method for borehole muon radiography orgeotomography, the method comprising: positioning a plurality ofborehole muon detectors in one or more boreholes, each borehole muondetector including a plurality of sealed drift tubes to form a bundle ofdrift tubes, each drift tube including a centrally located anode wiredisposed on a longitudinal axis, an inner surface which is coated with acathode coating, and a timer for measuring the time of pulses on theanode wire; generating a radial electric field between the anode wireand the cathode coating; and measuring a plurality of muon trajectories,wherein measuring the plurality of muon trajectories comprises, for eachpassing muon, determining X-Y measurements in a transverse plane of thepassing muon based on drift times of multiple drift tubes in the bundle.9. The method of claim 8, wherein determining the X-Y measurements inthe transverse plane of the passing muon based on drift times of themultiple drift tubes in the bundle comprises determining, for each ofthe multiple drift tubes, a radial position of the muon passing by theanode wire, and, based on the radial positions for at least three drifttubes, an azimuthal angle in the transverse plane based on fitting aline through isochrones of the at least three drift tubes.
 10. Themethod of claim 8 wherein measuring the plurality of muon trajectoriescomprises, for each passing muon, determining a zenith angle of the muontrajectory with respect to a longitudinal axis based on regression oflongitudinal positions of the passing muon measured in each of themultiple drift tubes.